Resonant low-energy electron scattering on short-range impurities in graphene 
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Resonant scattering of electrons with low energies (as compared to the bandwidth) on a single 
neutral short-range impurity in graphene is analyzed theoretically, taking into account the valley 
degeneracy. Resonances dramatically increase the scattering cross-section and introduce a strong 
energy dependence. Analysis of the tight-binding model shows that resonant scattering is typical for 
generic impurities as long as they are sufficiently strong (the potential is of the order of the electron 
bandwidth or higher). 



I. INTRODUCTION 

Electron transport in graphene is a subject of intense 
study, both theoretical and experimental, since the very 
discovery of this material in 2004. A In general, electron 
transport is determined by competition of different scat- 
tering mechanisms, both inelastic (e. g., electron- phonon) 
and elastic (static defects) . Elastic scattering is dominant 
at sufficiently low temperatures. 

Different kinds of crystal imperfections can cause elas- 
tic electron scattering in graphene: mesoscopic corru- 
gations of the graphene sheet (ripples) 2 - producing per- 
turbations smooth on the atomic scale, charged impu- 
rities producing long-range Coulomb fields, dislocations 
producing long-range strain fields, or short-range neu- 
tral impurities. While the first three types seem to be 
more important for the transport in the clean graphene, 
it is probably the fourth one that can be controlled. In 
Refs. y|4j graphene oxide was chemically reduced to nor- 
mal graphene; in Ref. [H hydrogen and hydroxil groups 
were deposited on the graphene sheet in a controlled and 
reversible manner; in all cases resistance changed by sev- 
eral orders of magnitude. Assuming that attachment of 
a chemical group to a carbon atom in graphene changes 
the hybridization of its electronic orbitals from sp 2 to 
sp 3 , one can view such group as a neutral short-range 
impurity. 

Short-range impurities in carbon nanotubes have 
been studied even before the single-layer graphene 
was obtained in the laboratory^^ 9 -^ Short-range im- 
purities have been shown to modify local electronic 
properties of graphene, such as the local density of 
state a 1:L i 12 ' 13 i 14 ' 15 i 16 ' 17 or local magnetic moment) 13 i 18 
and to induce Friedel oscillations in doped samples. 19 
The present work is dedicated to the problem of elec- 
tron scattering on a single short-range impurity whose 
size R is assumed to be of the order of the interatomic 
distance a (the C-C bond length), and the electron en- 
ergy e is assumed to be much smaller than the energy 
scale set by the potential, v/R (v is the electron veloc- 
ity at the Dirac point, so that v/R is of the order of the 
electronic bandwidth). The main focus is the case of a 
strong impurity so that the results for electron scattering 
obtained in the first Born approximatio n 20 ^ 1 ^^ 3 -^^ 
are not expected to be applicable. Instead, we are going 
to exploit the smallness \e\R/v <C 1. For particles with 



parabolic spectrum such low-energy scattering is char- 
acterized by a single constant of the dimensionality of 
length (the scattering length), determined from the so- 
lution of the Schrodinger equation at zero energy. In 
graphene, due to the degeneracy of the spectrum at the 
Dirac point, more than one length is needed to charac- 
terize a scatterer^ 

The common intuition is that for a strong enough scat- 
terer the typical value of scattering lengths I ~ R, yield- 
ing the cross-section a ~ pi 2 <C a (here p — \e\/v is the 
electron momentum counted from the Dirac point). The 
exception to this is the case of resonant scattering when 
the potential has a (quasi-)bound state with small en- 
ergy, then one of the scattering lengths becomes of the 
order of the size of this state. The main motivation for 
the present study is that for the Dirac spectrum the ex- 
ception becomes a rule: a vacancy (which can be viewed 
as the limit of an infinitely strong scatterer) introduces 
a bound state exactly at the Dirac point jii so one of the 
scattering lengths diverges. As a consequence, the scat- 
tering cross-section diverges as the electron energy ap- 
proaches the Dirac point. This divergence corresponds 
to that found in the scattering off vacancie o 13 ' 27 and in 
the unitary limit of potential impurities^ 

It is hard to introduce a real vacancy in graphene, how- 
ever, if the 7r-orbital of some carbon atom is very tightly 
bound, this atom acts effectively as a vacancy for the rest 
of the 7r-electrons in the crystal. In particular, impuri- 
ties introduced electrochemical!}* 3 ^^ are likely to act as 
strong scatterers, and thus are unlikely to be described 
by the first Born approximation. This expectation is sup- 
ported by density-functional theory (DFT) calculations 
for graphane^ 9 . (a hypothetical material obtained from 
graphene by attaching a hydrogen atom to each carbon) : 
the binding energy was obtained to be about 7 eV per 
hydrogen atom. For such scatterers the scattering length 
should be large, thus motivating the present study. 

The main framework of this study is the general scat- 
tering theory 3 ^ modified for the 2D Dirac equatio n 31 ' 32 i 33 
taking into account the valley degeneracy. In Sec. IHII we 
show that for a general short-range scatterer all the in- 
formation necessary to determine the cross-section up to 
corrections of the order (pR) 2 is encoded in a 4 x 4 ma- 
trix L and a constant ro (all having the dimensionality of 
length). The terms of the order (pR) 2 and higher cannot 
be studied using the Dirac equation, as the Dirac hamil- 
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tonian itself is the leading term in the expansion of a 
microscopic hamiltonian in the parameter pa (and R ~ a 
is always assumed here). The matrix L (i) can be ob- 
tained from the solution of the microscopic Schrodinger 
equation (e. g., an ab initio calculation) for electrons in 
the graphene crystal at zero energy and with an appro- 
priate asymptotics; (ii) is hermitian and invariant with 
respect to the time reversal, so it depends on 10 real pa- 
rameters; (iii) its four eigenvalues h, . . . , U play the role 
of the scattering lengths. Divergence of one or several 
of these eigenvalues signals the existence of a localized 
solution at zero energy. For the parameter r$ we have 
(i) r ~ R; (ii) the dependence of the scattering ampli- 
tude on rg is weak (logarithmic); (iii) the exact value 
of ro cannot be extracted from the zero-energy solutions 
only, wave functions at low but finite energies have to 
be considered in order to determine it. When L and ro 
are known, the low-energy scattering T-matrix is given 
by Eq. (|19[) . which covers both Born limit [T(e) is e- 
independent, the cross-section a cx e] and the unitary 
limit [T(e) cx l/(elne), a oc l/|eln 2 e|], as well as the 
crossover between them for an impurity of a large but 
finite strength. 

In Sec. lIVI wc consider two examples of impurities with 
special symmetries, and see how these symmetries man- 
ifest themselves in the scattering (i. e., how they restrict 
the form of the matrix L). The first example (the site- 
like impurity) is an impurity localized around one of the 
carbon atoms and preserving its C^ v symmetry (Cs v con- 
sists of three-fold rotations and reflections in three planes 
perpendicular to the crystal plane); it is natural to as- 
sume that this would be the case for a hydrogen atom 
bound to a carbon. The second kind (the bond-like im- 
purity) involves two neighboring carbon atoms and the 
bond between them and has the symmetry Ci v . This 
could be the case for an oxigen atom bound to two car- 
bon atoms. The two kinds of impurities, described above, 
are schematically shown in Fig. [1] Of course, a generic 
impurity is not going to have any symmetry at all. 

Since graphene crystal is symmetric, impurities of the 
same kind can occur in different locations and with dif- 
ferent orientations with equal probability if these can be 
related to each other by a crystal symmetry operation. 
For example, the site-like impurity, located on an A atom 
in Fig. [TJ can reside on a B atom with the same proba- 
bility; the bond-like impurity can have one of the three 
different orientations, rotated by 27r/3 with respect to 
each other. Although equivalent from the crystal sym- 
metry point of view, such impurities will have different 
L matrices. If one is not looking at effects of coherent 
scattering off several impurities, the cross-section can be 
averaged over such equally probable impurity configura- 
tions. This procedure is described in Sec. [V] 

In Sec. IVII we perform explicit calculations in the 
tight-binding model as an example of microscopic model 
(i. e., well-defined at short distances), using the stan- 
dard T-matrix approach on a lattice, adopted by many 
authors. 12 ' 13 i 14 ' 15 i 16 ' 17 We consider the two kinds of im- 
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FIG. 1: (Color online). The honeycomb lattice with two 
atoms (^4 and B) per unit cell and the elementary translation 
vectors ai and &2- The two inequivalent reflection planes a a 
and a' a are shown, (i) and (ii) are schematic representations 
of site-like and bond-like impurities with symmetries Cs v and 
C21), respectively. 



purities, mentioned above, modeled as a diagonal on-site 
potential for a site-like impurity, and a combination of 
a diagonal and an off-diagonal potentials for a bond-like 
impurity (the same model was adopted in Ref. |15T ). The 
scattering lengths are calculated as functions of the im- 
purity strengths, and ro is obtained to be 0.5 a. In agree- 
ment with the results of Ref. [HI, the divergence of the 
scattering length occurs at infinite impurity strength for 
a site- like impurity (corresponding to a zero-energy state 
bound to a vacancy) , and at finite values of the diagonal 
and off-diagonal strengths for a bond-like impurity. The 
resulting cross-sections as functions of the electron energy 
for different impurity strengths are shown in Figs. 12131 



II. FREE ELECTRONS IN GRAPHENE 

Graphene unit cell contains two atoms, labelled 
A and B (Fig. [TJ. Each of them has one 7r-orbital, so 
there are two electronic states for each point of the first 
Brillouin zone (the electron spin is not considered in the 
present work). The electronic energy e (measured from 
the Fermi level of the undoped graphene) vanishes at the 
two Dirac points K, K' with wave vectors ±K. Thus, 
there are exactly four electronic states with e = 0. An 
arbitrary linear combination of them is represented by a 
4-component column vector ip. Here we choose the fol- 
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TABLE I: Irreducible representations of the group Cev and 
their characters. 
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valley-diagonal matrices 
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A z 






Y^x-, Sy 


SyA 2 , Yj x A. z 




valley-off-diagonal matrices 




matrix 


S z A x 


Xj Z Ay 


A x 


Ay 




Sx-Ay, SyAy 



TABLE II: Classification of 4 x 4 hermitian matrices by irre- 
ducible representations of the Cev group. Matrices joined by 
braces transform through each other under translations. 



lowing arrangement of the wave function components in 
the column^ 



IpAK 
IpBK 
IpBK' 
-1pAK> 



(1) 



Other definitions of the column vector are possible, but 
they are inconvenient for the present problem; for the 
discussion see Appendix lAl and Ref. [H. Being interested 
in low-energy states, we focus on states in the vicinities 
of the Dirac points. The wave vectors of these states can 
be written as k = ±K + p, where pa C 1 (here a « 1.42 A 
is the C-C bond length). Equivalently, states near the 
Dirac points are obtained by including a smooth position 
dependence ^(r), r = {x,y). 

The basis in the space of 4 x 4 hermitian matrices 
is formed by 16 generators of the SU(A) group. They 
can be represented as products of two mutually commut 



the £ matrices act on the x variables, while the A ma- 
trices act on the y variables. The basis in the valley 
subspace is denoted by 4>k , 4>K' for future convenience. 

The matrices Si,Aj, i,j — x,y,z, and their products 
have definite transformation properties under the crystal 
group. The irreducible representations of C§ v group are 
listed in Table U (C n denotes the rotation by 2-7r/n, a a .b.c 
are the three reflections which swap the K and K 1 points, 
a 'abc are t ne three reflections which swap A and B atoms, 
the fixed point of these operations being the center of the 
hexagon, see Fig.[T|). The correspondence between them 
and the matrices is given in Table HI1 35,36 

The unitary matrices corresponding to the symmetry 
operations of the crystal can also be written in terms 
of £ and A matrices, independently of the representa- 
tion usedi^ Specifically, Uc 3 = e^ 2Tr ^ 3 ^ Ez is the ma- 
trix of the C3 rotation, Uc 2 — S Z A X - of the C2 rota- 
tion, U a ' a — T, X A Z - of the a' a reflection, U„ a = T, y A y 
- of the a a reflection. The two elementary transla- 
tions by the vectors aj., a 2 act on the wave function as 
t ai 2 : tp(r) i-> e T( - 27r ^ 3 ^ Az Tp(r — ai j2 ). The time reversal 
operation is defined as ip <— > Utip* , where the unitary time 
reversal matrix assumes the convenient form Ut = S y A y 
in representation ([!}. 

Taking the leading-order term in the expansion of the 
band hamiltonian in the powers pa < 1, we describe the 
electrons by the Dirac hamiltonian 



H 



■!«S • V, 



(3) 



where S = (T, x ,H y ) is a two-dimensional vector, and 
v w 10 s cm/s is the electron velocity. The eigenstates of 
the Dirac hamiltonian with a definite value of momentum 
are plane waves: 



J, 



(o) 

ifSK, 



1 

71 



*P r 7 /,(0) 



v 2 



(4a) 



(4b) 



with the energy e psK = svp. The index s = ±1 distin- 
guishes between the conduction and the valence band. 



ing algebras of Pauli matrices denoted by E x ,E y ,H z and The unitary matrices W Vv = e ffiWV s *<<W 2 , wne re 



A x , A y , A 2 ) 23 ' 26 which fixes their algebraic relations. We 
denote the unit 4x4 matrix by 1, and sometimes So or 
Ao to make the formulas compact. In representation (|TJ) 
Tii are the Pauli matrices acting within upper and lower 
2-blocks (the sublattice subspace), while Aj are the Pauli 
matrices acting in the "external" subspace of the 2-blocks 
(the valley subspace) . For 4-column vectors which can be 
represented as a direct product 



X2 Vi 



x-iyi 




X\ 




Vi 




X\ 


X1V2 




X2 




V2 




X2 



if p = &Yct&n(p y /p x ) is the polar angle of the vector p, 
diagonalizc the Dirac hamiltonian in the momentum rep- 
resentation: vp-H — Wl vpT, z W iPp (note that W v= 2-n = 
—W v —o). The index k = K,K' labels the valleys. As 
hamiltonian ([3]) does not contain A matrices, the valley 
subspace is degenerate, so any other basis can be chosen. 

Besides plane waves, we will need the wave functions 
of states with a definite half-integer value of the "total 
angular momentum" j z — —i(d/dip) + (l/2)'S z which also 
commutes with the Dirac hamiltonian ([3]): 



(yi^K + V2<t>K>), 



(2) 



sJ 32 _ 1/2 (pr)e l(3 ^ 1/2)(v+T/2) 



J 



J»+l/2 



(pr)e 



i(j z +l/2)( V +n/2) 



(5) 
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where J m are Bessel functions of the first kind. If one replaced by the Bessel functions of the second kind Y T 
relaxes the condition of regularity of the wave function or Hankel functions h£ 2) = J m ± iY m . 
at r = 0, the Bessel functions of the first kind J m can be 



III. RESONANT SCATTERING ON A SHORT-RANGE POTENTIAL 



General definitions 



Here we use the standard expansion in partial waves^ modified for the Dirac equation analogously to Refs. l3lll32ll33l . 
For a potential V(r) that falls off rapidly at distances r > R, the electron motion at distances r 3> R can be considered 
free. The general scattering solution corresponding to the energy e = svp can be written as 



Vw(r) = 

rn=—o 

+ E 4 



, — i(m+l/2)v3 p 



V2 



sJ TO (pr)e im ^ +7r / 2 ) 
J m+1 (pr)e^ m+1 ^ +7! ^ 



v/2 



rC 1 )/ — \ p im(<?-Hr/2) 



sH^(pr)e l 



i(m+l)( ¥ j+7r/2) 



1 



•7^+1/2 (P)<^- 



(6) 



The first sum represents just the incident plane wave (J1J). The second sum in Eq. © represents the outgoing scattered 
wave. Since a short-range potential can change arbitrarily the structure of the state in the valley subspace, an arbitrary 
2x2 matrix -^+1/2 (p) ^ s introduced (the factor 1/2 in front of T is introduced for convenience)^ The 2x2 matrix T 
can be represented as a linear combination of A matrices. 

Using the asymptotic behaviour of Hm (pr) at pr 3> \m 2 — 1/4| for integer m: 



'2/tt 
pr 



} ipr—im7r/2 — i7r/4: 



(7) 



and the standard definitions of the scattering amplitude f s (<p, Vp'iP) an( i the scattering matrix S s ((p,ip';p) in two 
dimensions j^l where the existence of the degenerate valley subspace is taken into account: 



Vw( r ) 



^ Pr < } S K 

3 — ipr-\-i7T / 4 



3 ipr+27r/4 



y/2irpr 



27r^ r -^p-7r)4 ) + 



gipr— i-n/A 



yJ2npr 



S s (<p T ,<p p ;p) ip$ aK , 



(8a) 
(8b) 



we relate them to the matrix .T^+i/aCP)' calculating the probability current j(r) = ip'(r)v'Eip(r), we obtain the 
differential cross-section dcr SKK i(ip, <p';p) for an incident particle with momentum p = (jpeosep' ,p sin ip') in the valley k' 
to be scattered into the valley re in the direction n = (cos ip, sin ip): 



fs(<p, <p';p) 



T s (ip,ip';p) 



i^fhrp 

OO 

7 

da SKK >(tp,ip';p) 



m— — oo 



dip 



\4>kfs{<P, V>'\P) <t>K' 



(9a) 
(9b) 
(9c) 

(9d) 



Note that VSf27r,«« = the above definitions imply < <p T — tp p < 2n. For — 2tt < ip r — ip p < we have to set 

S a (ip,tp';p) — —S s (tp + 2tt, ip'; p). The scattering matrix satisfies the unitarity and reciprocity conditions (the latter 
assuming the symmetry of the scattering potential with respect to the time reversal): 



Z7T 

f dip 

o 

S s ((p,(f';p) = k y Sj(tp' +TT,ip + Tr;p)A y , 



f2), 



(10a) 
(10b) 



5 



where S T denotes the transpose of the matrix S. The scattering amplitude can be related to the T-matrix T(p, p'; e) 
on the mass shell, |p| = |p'| = \e\/v. Starting from the exact expression for the scattered wave function, 

^ K (r) = e^<L + J |^ e'P' r G(p', svp) T(p', p; svp) G( P , e) ee /^^ P '^ +)2 , (11) 

and taking its ^> 1 asymptotics, we arrive at [n = (cos if, sin (p) and n' = (cost/?', sin <£>') are unit vectors]: 



WpT(pn,pn , ;M;p) Wj, = 



it) 1 
?2 



1 + s 
1- s 



sF s {tp,ip')p). 



(12) 



B. Scattering lengths 



At 



R the potential mixes different terms in Eq. 



The asymptotic behaviour of the Bessel and Hankel 



functions at pr <C \/\m\ + 1 for integer m is given by 

(signm)"' /pr\l m l 

Jm{pr) 

(i) 



H 



(1)/ 



\m\\ \ 2 

(signm) m |m|! / 2 
i7r|m| \pr 



\m\ 



l/«') = l + |(ln^+ ). 



(13a) 

(13b) 
(13c) 



where 7 = 0.5772 ... is the Euler-Mascheroni constant. At this stage we make an assumption that all terms constituting 
the scattered wave in Eq. ([6]) should be of the same order at r ~ R, provided that their coupling is allowed by the 
symmetry of the scattering potential. Depending on the symmetry of the scattering potential, we have to consider 
two cases, (i) The potential is isotropic and S z -conserving: V{v) — V(r), Y> z V(y)Y> z = V(r) (this case was analyzed 
in Refs. l3lll32ll33h . Then j z = —i(d/dtp) + (1/2)S Z is conserved, so the terms in Eq. © with different values of m are 
decoupled (each term with a given m corresponds to j z =m+l/2). Matching the terms gives Tj a ~ (pi?) 2 '^'. (ii) The 
potential V(r) is generic, so it mixes all states with different values of j z . Matching the terms gives J r j z ~ (pR)^ z ^ +1 ^ 2 . 
Thus, F±i/2 ar e the most important terms in both cases; moreover, as typically R ~ a, considering the terms with 
\j z \ > 1/2 would require going beyond the Dirac hamiltonian, since the Dirac hamiltonian itself is the leading term 
in the expansion in pa. 

At r <C 1 jp we can neglect the energy in the Dirac equation, which becomes 



[-ivH ■ V + V(r)]i/}(r) = 0. 
At r ^> R this equation admits solutions of the form 



(14) 



(15) 



which determine the asymptotics of different angular harmonics of the four linearly independent zero-energy solutions. 
Since non-zero angular harmonics of the incident wave in Eq. © vanish at p — > 0, we are interested in the solutions 
of Eq. (fT4| whose asymptotics can be writen as 



1 




e iip 
ir 



1 





e itp 





1 


ir 


1 



m> 1 



o 



L220K 

m>l 



(16a) 
(16b) 



where each Lij is a 2 x 2 matrix in the valley subspace, which has to be determined from the solution of the Schrodinger 
equation at short distances (the factor 1/i is introduced for convenience). Let us associate the indices i,j = 1,2 of 
the 2x2 matrices with the matrix sturcture in the S-subspace, thus combining the four 2x2 matrices into 
a single 4x4 matrix L. Then the asymptotic behavior of any solution of Eq. (TTJJ at r — > 00 can be written as: 



1 



^(r) = ^> + — (n • S)i^> 
ir 



(0) 



m> 1 



(17) 
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where is an arbitrary 4-column. This equation could also be viewed as the boundary condition on the angular 
harmonics of the scattering solution ^ at r — > 0; being formed at short distances r ~ R, this boundary condition 

should not depend on e for |e| <C v/R. However, due to the logarithmic divergence of Hankel function Hq (pr), 
matching of wave functions should be performed at some r = ro ~ R. Note that in contrast to the two-dimensional 
Schrodinger equation, the value of the constant ro cannot be determined from the zero-energy solution since the 
logarithmic function is not a solution of the Dirac equation at zero energy. 

Comparing expressions l|16j) to the scattering solution ([6|) and using asymptotic expressions (| 13[> . we obtain the 
general possible form of F±i/ 2 (p)'- 



— -?+i/2(p) = sL n 

np T ' 



-^/ 2 + iffW(pro)^ 1/2 (p) 



np 



■T s _ 1/a (p) = i 2 i 



-itpp/2 



4Vo)ni/2(p) 



l-iia 



; W2 + l if m (pr . o)j r :j 2(p) 



(18a) 
(18b) 



Solving these equations and comparing the result to 
Eq. I|12p , we obtain the general low-energy T-matrix4£ 



IV. IMPURITIES WITH SPECIAL 
SYMMETRIES 



T(e) 



In 



2v 



?,7T 

y 



2ttoL, (19) 



which determines the differential cross-section: 



dip 



2nv 2 



(^jT(svp)^X' ■ ( 20 ) 



The matrix L satisfies (i) L = to ensure the unitarity 
of the scattering matrix IjlOap , and (ii) L = UtL T [lJ as a 
consequence of the reciprocity condition (|10b[) . Thus, 
it has four orthogonal eigenvectors ipi~. Lipi — Uipi, 
i = 1,...,4, and the four eigenvalues U play the role 
of the scattering lengths. The angular dependence of 
dcr SKK '(tp,(p';p)/dcp is given by the sum of an isotropic 
term and terms cx e ±lv , e ±l<p ' , and e ±ltp±lv ' (with all 
four combinations of the signs). The total out-scattering 
cross-section (i. e., integrated over ip and summed over v) 
takes the simple form: 



4 



spln(eipr /2)] + {iTp/2f 



(21) 

Typically, one assumes that for a strong potential all 
scattering lengths U ~ R. However, an explicit calcula- 
tion for a point defect in the tight-binding model, per- 
formed in Sec. IVI1 shows that one of the lengths U (let it 
be li for definiteness) can become arbitrarily large. The 
case l\ —* oo corresponds to the existence of a localized 
solution ip(r) ~ (n- S)-0i/f* at zero energyAi In this case 
the cross-section diverges at p — > 0. This divergence cor- 
responds to a similar divergence in the i mag inary part of 
the electron self-energy found in Rcfs. [l3ll28c a similar di- 
vergence in the cross-section was found in Ref. [I?]. Note 
that even in the case of resonant scattering the scaling 
Tj z ~ {pR)^ z ^T±i/2 holds: indeed, at l\ — * oo the co- 
efficients at the l/r m terms in the wave function of the 
localized state should scale as R m , as there is no other 
length scale in the problem. 



Let us consider two particular kinds of impurities, 
shown in Fig. [T] 

(i) A site-like impurity with the symmetry C% v whose 
fixed point is located on one of the atoms (let us assume it 
to be an A atom). Thus, the matrix L should be invariant 
under the reflection a a and the rotation C3 = 6*3^1 (we 
remind that the rotation C3 is around the center of the 
hexagon). The conditions L = U^ a LU aa , L = U^,,LUc> 3 
together with the time-reversal symmetry restrict the 
matrix L to 

L = La 1 1 + L B2 A 2 S 2 + L E (A X E X - A y Z y ). (22) 

The eigenvalues and eigenvectors of this family of matri- 
ces are 

h.2 = La 1 + Lb 2 ± 2Le, h,4 = La ± — Ls 2 (23a) 



V2 



1 



±1 



, ft 



3,4 



V2 





1 

±1 





(23b) 



The T-matrix has the same transformation properties as 
the matrix L, so it can be written in the same form (|22j) 
with the substitution L — > T. Then, according to 
Eq. (|2T)1) the differential intravalley and intervalley cross- 
sections can be written as 



da 



KK 



P 



2ttv 2 



dip 

d<J K 'K V 



dtp 



2nv- 



Ta 1 cos 



tp-tp 



iTs 2 sin 



tp-ip 



(24a) 
(24b) 



(ii) A bond-like impurity with the symmetry C 2 ^ whose 
fixed point is located at the center of a bond (let us as- 
sume it to be a bond connecting the two atoms within the 
same unit cell). Then, the matrix L should be invariant 
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under the reflection a' a and the rotation C 2 = Citg^t^, 
which fixes 

L = Lax^ + Le 2 ^- z ^"x + Lax^-x^z + LE 2 A y '£,y, 



h,2 — (La! + L E2 ) ± (L E . 2 + L Al ), 
h,A = (Lai — Le 2 ) ± (Le 2 — Lax), 



1pl,2 = 



1 



±1 
1 
1 



^3,4 



1 

±1 

Tl 
1 



(25a) 
(25b) 
(25c) 

(25d) 



The differential intravalley and intervalley cross-sections 
are 



dip 

ddK'K 

dip 



P 



2ttv 2 



2ttv 2 



Ta x cos 



sTe 2 cos 



<p + <p' 



±Ax sin — h si e 2 sin — - — 



(26a) 

2 



(26b) 



If the location of the impurity is different from what we 
have assumed, its L matrix can be obtained by applying 
the corresponding symmetry operation. For example, in 
case (i) the matrix for an impurity located onaB atom is 
obtained by a C2 rotation: L — ► A X T, Z LA X Y. Z . Obviously, 
these two locations can occur with equal probability, so 
one could average over them. This procedure is described 
in the next section. 



V. AVERAGING OVER THE IMPURITIES 

As discussed in the end of the precedeing section, the 
presence of defects characterized by a certain T-matrix 
T(p, p', e) implies the presence of the same (on average) 
number of defects of the same type, placed in differ- 
ent locations with different orientations, and thus hav- 
ing different T-matrices, but equivalent with respect to 
the symmetry of the crystal. If one studies effects which 
do not involve coherent scattering on several impurities, 
it is sufficient to average any observable 0[T(p, p', e)], 
calculated for a single impurity, according to 



O 



3\Cev\ 



]T (o[u R T(Rp,Rp',e)U 



Here \Cq v \ — 12 is the number of elements in the 
Cq v group, R are the operations from the group, and Ur 
are their 4x4 matrices in the -0-representation. The aver- 
aging is performed also over the elementary translations 



with the matrices Ut. 



It is convenient to consider the group Cg v - the direct 
product of the point group C'q v and the 3-cyclic group 
represented by the matrices 1, e ±( 27 ™/ 3 ) A » _ Then Eq. (j2"?| 
describes simply the average over the group Cq v . 

Let us apply this procedure to the differential cross- 
section. We write the averaged Eq. (|2"0|) as 



da 

dip 



P 



2irv 2 



E 



Tr{ UU V-V* ) UrT* U r ) U R T} 



\CL\ 



(28) 

where we abbreviated ip = ip$] K , ip' = ipffl SK >, and 
T = T(svp). Thus, equivalently, we can calculate the 

average (ip'ip'^) ® (ipip^). In the matrices ipip^ and ip'ip'' 
we separate the components corresponding to different 
irreducible representations of Cq v : 

= \ J2 (^AW) AiEj = 

i,j=0,x,y,z 

1 



(1± A z )(l + sn- £), 



(29) 



where the plus (minus) sign should be taken for <p K = tpK 
(<p K = 4>k'), and n = (cos ip, sin ip). Let us label the 
matrices AjEj belonging to an irreducible representa- 
tion r of the dimensionality d r as (AS) J, where the in- 
dex £ — 1, . . . , rf r labels the matrices within the repre- 
sentation. Then in each representation we can define the 
d r x d r matrices {U R )ui as 



4(A£)^=X>^)«,(A£)^ 



(30) 



'=1 



From the orthogonality relation for the representation 



matrice: 



,38 



1 ^ — ^ ' ^rr' 

TFT - ! ( U R)iil2( U R)i3iA = ~T faits^lt, ( 31 ) 



O 



Ul a Ur T(Rp, Rp', e) U R Ut ai J . (27) we obtain the general expression: 
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da SKK >{<p,ip';p) 
dip 



P 



Tr{AoTi{svp)A T{svp)} + ^cos(<p-<fi') ^ Tr {A ^ j T^(svp)A E j T(svp)} 



32ttv 2 
+ {25 KK , - 1) [A -> A 



(32) 



r 



(here we denoted the 4x4 unit matrix 1 by Ao). 

To conclude this section, we note that the averaging 
procedure described above is equivalent to averaging the 
cross-section da{(p, if') j dip over (tp + tp')/2 keeping <p — (p' 
fixed (due to the symmetry of the crystal with respect to 
C3 rotations), and subsequent averaging over the sign of 
tp — ip' (due to the symmetry with respect to reflections). 



VI. IMPURITIES IN THE TIGHT-BINDING 
MODEL 

A. Green's function 

Let us consider the tight-binding model with nearest- 
neighbor coupling as an exactly solvable example of 
a microscopic model (i. e., well-defined at short dis- 
tances). The only parameter of the clean hamiltonian 
is the nearest-neighbor matrix element which we write 
as —2v/(3a), thus expressing it in terms of the electron 
velocity at the Dirac point. It is convenient to work with 
a 2-component wave function {^^(r„), 5 , s(r n )} (corre- 
sponding to the two atoms in the unit cell), where the 
position of the unit cell r„ = nxa.% + n 2 &2 is labelled 
by two integers ni,ri2. The tight-binding hamiltonian 
H.o(r n — r n /) is a 2 x 2 matrix in the sublattice space. 

The scattering problem in the tight-binding model 
with a few-site potential 14 is conveniently solved using 
Lippmann-Schwinger equation: 



(33) 



where is the incident wave, ^> is the sought wave 
function, and Q(e) = (e — Ho)~ 1 is the Green's function, 
explicitly given by 



d 2 k e 4k ( r "- r - 



A BZ e 2 -|tk| 2 



e 



(34) 



3a v 



-ikai 



ika 2 



A B z = 



(2irf 



(27T) 



27a 2 /2 



The large-distance behavior of Q(r, e) is determined by 
the singularities of the denominator, i. e., vicinities of the 
Dirac points ±K, k = ±K+p, where we can approximate 



1 



-i(±K+p)ai 



-i(±K+p)a 2 



3« 



(TPx + iPy)- (35) 



Focusing at |e| <C v/a, we obtain for r ^S> a (c. c. stands 
for the complex conjugate): 



GAA{r,e) 
Gba{y, e) 
Gba(t,0) 



Aickl /iKr 

Nuclei 



(V K 



Aiv 2 

^uc e 



V 



2-Kir 



e- Kr ) sH^(\e\r/v), 

(36a) 

' HP(\e\r/v), 
(36b) 

(36c) 



r -\-i(p-\-iir f 2 



-iiCr—iip 



We also need the Green's function at coinciding points: 
SAA(0,e) = --5- (in— - 7+ -)+0(e 3 ), 



|e|r 



G B A(0,e) 



2v 



0( e 2 ). 



(37a) 
(37b) 



The value of tq in Eq. (|37a|) is determined by the in- 
tegration over the whole first Brillouin zone; numerical 
integration gives e 7 ro = a within the numerical precision. 
The leading term in Eq. (|37b[) can be easily obtained in 
the coordinate representation using the fact that H^ 1 , 
just like H.q, is invariant under C3 rotations around each 
carbon atom. 



B. One-site impurity 

Let us add the on-site potential Uq different from zero 
only on the A atom of the n\ — n 2 — unit cell. The 
limit Uo — > 00 is equivalent to imposing the boundary 
condition ^a(0) = and thus describes a vacancy. 

First of all, we note that the two plane wave states with 



*A(r„) = 0, *s(r n ) 



3 ±iKr„ 



remain the zero-energy 



eigenstates of the hamiltonian even in the presence of 
the potential. In the representation ([I]) these two states 
are represented by the 4-columns 



'K' 



Comparing them to the asymptotic forms (fTB)) . we see 
that L\\4>k> = L 2 i(j>K' = 0, L 12 <f>K = L 22 cj) K = 0. 



The other two zero-energy solutions correspond to the 
incident wave on the A-sublattice, f ^(r n ) = e ±jKr ™, 

(r n ) = 0. In the representation ^ these two states 
are represented by the 4-columns 
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-1 





-1 



As the potential Uo is localized on one atom, and 
G(e = 0) is off-diagonal in the sublattices, the Lippmann- 
Schwinger equation is straightforwardly solved to give the 
wave function: 



* A (r„) = e 



±iKr„ 



<Mr„) = C/o£sA(r„,e = 0). (38) 



Using Eq. (|36cp . we obtain 

L = A UC U 1 + A Z T, Z + A y Z y - A X E X 



(39) 



U ( = 10eV 




-10 12 

electron energy (eV) 

FIG. 2: The out-scattering cross-section (coinciding with the 
transport cross-section) for a one-site impurity of the strength 
Uo = 1,5, 10 eV (dotted, dashed, and solid curve, respec- 
tively) as a function of the electron energy (in eV), as given 
by Eq. (|42|> . The curve for Uo = leV is multiplied by a factor 
of 10. The parameters of the model are v = 10 8 cm/s = 6.58A, 
a = 1.42 A. 



in agreement with Eq. ((22)) . The eigenvalues of this ma- 
trix are easily found to be ij = A uc U /(nv), h — h — 
I4 = 0. At Uo — > 00 the scattering length l\ diverges. In 
this case the amplitude of the incident wave can be set to 
zero, and \ffl(r n ) oc Gba{?, e — 0) is the wave function 
of the state, localized on the vacancy. 

At e 7^ the Lippmann-Schwinger equation is solved 
self-consistently for ^^(0) to give the wave functions: 



C. Two-site impurity 



Let us add a potential which mixes the two sites in the 
m = ri2 — unit cell: 



*k s (r„) = 



i* k /2 



+ 



,(40) 



Uq 1 - GAA(0,e) \GBA(r n ,e 
where e** k = ik/|tk|< This corresponds to the T-matrix 

1 + A Z T, Z + AyEy - A X T, X 



T(e) = 



Uq 1 — Gaa(0, e) ' 



(41) 

Using Eq. ([37a]) . we arrive at Eq. (13]). 

To calculate the average cross-section, we note that in 
Eq. (]32[) only the first term survives, so the scattering is 
isotropic in space and completely mixes the valleys. The 
total out-scattering cross-section and the transport cross- 
section, averaged over the impurity positions, coincide 
and are given by 



a out (e) 



w 2 v\e\/2 



U A U 



hi 



2v 



e 7 r |e| 



(42) 



U = 



Uo Ui 
Ui Uo 



(43) 



We have chosen U\ to be real in order to preserve the 
A «-> B symmetry. The self-consistent solution of the 
Lippmann-Schwinger equation gives 



Wtf(O) = [U- 1 -G(0,e)] 1 VW(0) = T(e)¥<°>(0). 

Comparing the resulting wave function with Eq. (|17[) , we 
obtain 



2ltV 



1 + A z s 1 - A, 

no) + £„T(0)E„ — ^ 



/ n„ A T + iA„ „ , A T — iA„ 

iT(0)E a — y - + iZyT(0) — v - 



(45) 



This cross-section is plotted in Fig. [5] as a function of e 
for several values of Uo = 1, 5, 10 eV. 



For the potential of the form (|4"3")l the scattering lengths 
are obtained as [the matrix form of L is in agreement 
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-10 12 
electron energy (eV) 

FIG. 3: The out-scattering cross-section for a two-site impu- 
rity with Uq = 5 eV, U\ = (dashed curve), and Uo = 0, 
Ui = 5 eV (solid curve) as a function of the electron energy 
(in eV), as obtained from Eq. (|47a[) . The parameters of the 
model are v = 10 8 cm/s = 6.58 A, a = 1.42 A. 



with Eq. {25aJ] 
T(0) 



(£/ 2 - u?)u Q 



+ 



(t/ 2 



1 



[U 1 + (a/2v)(Ul-Ui)}- 
-U?)(U 1 + (a/2v)(U 2 -U 2 )) 



U 2 
1 o 



L = 



h, 



= %t 

2ttv 



-[Ut 



u ±u 1 



(a/2v)(U 2 -UfW 
A T 

^UC 



2-kv 



(A z E a 



ttw 1 - {a/2v){Ux ± E/ )' 



3,4 



(46a) 
■A a S ), 

(46b) 
0. (46c) 



The scattering lengths diverge when U\ = 2v/a ± Uq, in 
agreement with the results of Ref. [HI. 

Calculation of the T-matrix from Eq. (fTO)) and its sub- 
stitution into Eq. ([52")) gives the differential intravalley 
and intervalley cross-section (we set ip' = 0, as it depends 
only on <p — ip'): 



dcr KK > _ nv\e\ 
dip ~ 8 

*i,a(e) 



M 2 + M 2 ± 
i 



|*1±*2| 



■ cos </? 



u/ii,2 + e]n[2t!/(eTro| 



me/2 



(47a) 
(47b) 



In Eq. (|47a[) the upper and lower sign is taken for the 
intravalley (k = k') and intervalley (k ^ k') scatter- 
ing, respectively. The total out-scattering cross-section 
is plotted in Fig. [3] for the two cases of Uq = 5eV, U\ = 0, 
and U = 0, U x = 5 eV. 



VII. CONCLUSIONS 

In this paper we have studied scattering of low-energy 
electrons on a single neutral short-range impurity in 



graphene within the framework of the 2D Dirac equation 
taking into account valley degeneracy. We have shown 
that for a general short-range scatterer the most impor- 
tant information needed to determine the cross-section is 
encoded in a 4 x 4 matrix L, whose eigenvalues are the 
scattering lengths. Divergence of one or several scatter- 
ing lengths occurs whenever the impurity has bound elec- 
tronic states exactly at zero energy, which is accompanied 
by the singular behavior of the scattering cross-section as 
a function of the electronic energy. Quasi-bound states 
manifest themselves as resonances at finite energies, their 
width determined by the energy itself. 

The matrix L can be obtained from the solution of 
a microscopic model for the impurity in graphene; only 
the zero-energy states need to be considered for this. As 
an example of a microscopic model we take the tight- 
binding model and calculate the scattering lengths for 
the diagonal one-site impurity potential and the two-site 
potential having both diagonal and off-diagonal compo- 
nents. We obtain that one of the scattering lengths in- 
deed becomes much larger than the interatomic spac- 
ing for generic strong impurities (i. e., when impurity 
strength is of the order of the electronic bandwidth). 
This results in (i) a dramatic increase of the scattering 
cross-section and (ii) its strong energy dependence. 
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APPENDIX A: DIFFERENT 
REPRESENTATIONS OF THE STATE VECTOR 

If one adopts a representation, different from (TTJ) (let 
us denote the corresponding 4-columns by ?/>), the ma- 
trices Sj, Aj, i,j = x, y, z, instead of being simple Pauli 
matrices become some 4x4 matrices. Their algebraic 
relations and symmetry properties, listed in Table HT1 re- 
main the same. In fact, the convenient way to define 
these matrices for an arbitrary representation is to spec- 
ify the irreducible representation of the Cq v group, the 
point group of graphene, according to which they trans- 
form. This is sufficient to fix their algebraic relations. 34 
For example, the isospin matrices are defined as 

the matrices, diagonal in the K, K' subspace, and trans- 
forming according to the E\ representation of C§ v . 

As we have defined mutually commuting matrices Si 
and Aj, i,j = x,y,z, and have written the free electron 
hamiltonian in terms of the S matrices only, we must 
separate the degenerate valley subspace, defined as that 
invariant under the action of the S matrices. Since the 
basis vectors of representation {l} already have the neces- 
sary structure of the direct product, this representation 
is preferred for dealing with scattering problems. One 
can pass to it by choosing four basis vectors tpi, . . . , ip&, 
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defined as eigenvectors of £ z and A z : 

E^l = A z -0! = -0!, 
-£ z "02 = A z 2 = i>2, 

(Al) 

£^3 = -A z -03 = "03, 
-£ z 04 = -A z 04 = 04- 

Their relative phases are fixed by the requirement that 
the matrices £j act as the Pauli matrices in the subspaces 
{01,02} and {ip3,ipi}, and A* act as the Pauli matri- 
ces in the subspaces {ipi,ipa} and {02, Thus, vec- 
tors 0>i , . . . , 04 can be identified with the basis columns 
[10 0] T , [010 0] T , [0 010] T , [0 01] T in representa- 
tion JT]). up to an overall phase. The two representations 



are related by a unitary matrix U: tp = Uip. 

The overall phase of the matrix U is fixed by requir- 
ing the proper form of the time reversal matrix. In the 
-0-representation the unitary time reversal matrix U t can 
be different from H y K y . Indeed, the matrices in the two 
representations are related by U t = UU t U T (while Ej 
and Aj are transformed by applying U and J7 T ), which is 
sensitive to the overall phase of U. However, the prop- 
erties U t *U t = 1, and UtEfuJ = -E 4 , U t A*Ut = -A*, 
i = x,y, z, do not depend on the representation. Apply- 
ing these relations in the newly constructed representa- 
tion, we obtain Ut = E^Aj, up to a phase; this phase is 
nullified by the appropriate choice of the phase of U. 
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